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Abstract: VERITAS is an array of Imaging Atmospheric Cherenkov Telescopes designed for very 
high energy gamma ray (E>100GeV) observations of astrophysical sources. The experiment began its 
scientific observation program in the 2006/2007 observing season. We describe here the analysis chain 
for reducing the data, reconstructing the direction and energy of incident gamma rays and the rejection of 
background cosmic rays. 



Introduction 

The VERITAS instrument, described in more de- 
tail in [9], consists of an array of 4 telescopes each 
equipped with 499 pixel cameras [10] read out by 
500MS/S flash- ADCs (fADC) [5]. The traces are 
recorded into a custom data format before being 
archived for distribution to the collaboration and 
further analysis. Diagnostic plots - such as indi- 
vidual telescope and array trigger rates, deadtime 
stabilitity, atmospheric clarity, etc - are generated 
from these files and database stored information. 
These plots are made available through a web page 
for data integrity checks prior to an analysis be- 
ing attempted. A separate web page listing plots 
showing the long term history of the camera pixels 
(trigger participation, voltage information, etc) in 
order to trace camera health is also available. 

A number of packages have been developed for the 
analysis of VERITAS data. Having multiple, in- 
dependent packages available guards against any 
systematic errors being introduced into the analy- 
sis chain. Every analysis of a putative source must 
then have an independent secondary analysis in or- 
der to confirm any results. This paper describes the 
chain appropriate to the eventdisplay [6] and the 
VEGAS [11] packages. There are others, such as 
GrISU [2] and the quicklook package (used for on- 
line analysis when an observation is in progress) 
which show consistent results to the ones given 
here. 



Data Analysis 

The data analysis chain has three distinct elements: 
firstly the pixel data is calibrated to get accurate 
charge information; after calibration the pixel data 
is parameterised to give telescope level images, 
which in turn can be combined to reconstruct the 
geometrical properties of the shower; finally an es- 
timate of the background level of counts is made in 
order to determine if there are any statistically sig- 
nificant excesses (i.e. sources) in the field of view 
(fov). Once a source has been identified a data 
analysis can be further extended to look for vari- 
ability in emission with a light curve and a spec- 
trum calculated to try and distinguish between dif- 
ferent models of the emission process; or if there is 
no significant signal an upper limit can be placed 
on the source emission. Specific details of the 
higher level analysis process are left to the papers 
describing the analysis of specific objects. 

Calibration 

The technical details of obtaining calibration infor- 
mation are covered in [4] . Laser runs are used to 
calculate the relative gain between the pixels and 
the timing offsets to counteract any time spread 
due to path length differences (e.g. differing cable 
lengths, etc) for each channel. Pedestal events are 
injected at 1 Hz during an observation run to give 
an estimate of both the expected voltage offset for 
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each fADC trace and an estimate of the night sky 
background (NSB) and electronics noise from the 
variance of that value. VEGAS allows the pedestal 
variance to be calculated as a function of time into 
the run, with a default window of 3 minutes being 
the average time that a star will take to pass through 
a pixel's fov. Pixels with values outside of an ex- 
pected range for the calibration parameters are then 
flagged as bad and removed from further analysis. 
The charge values for these problem tubes are set 
to and typically this will apply to < 5% of chan- 
nels across the array (depending on the number of 
bright stars in the fov). 

The amount of charge in a good trace is then calcu- 
lated by summing the samples for a given window 
size. In VEGAS this is defaulted to a 14 ns integra- 
tion window that starts from a position calculated 
by averaging the pulse arrival time (defined as the 
half-height point on the pulse leading edge) of a 
representative number of events. In eventdisplay a 
two pass scheme is used. The first pass has a wide 
20 ns window that begins at a fixed position in or- 
der to calculate the integrated charge and the pulse 
arrival time. The second pass tightens the integra- 
tion window to 10 ns in order to increase the signal 
to noise ratio, placing the start of the integration 
according to the calculated time gradient across the 
image after it has been parameterised, as discussed 
in the next section. 

Parameterisation and reconstruction 

In the process of cleaning an image for parameter- 
isation, pixels which produce an integrated charge 
greater than 5 times their pedestal standard devi- 
ation (picture pixel) and any pixel that is adja- 
cent to these higher threshold pixels and having 
2.5 times their own standard deviation (boundary 
pixel) are automatically assigned to the image. Iso- 
lated picture pixels, i.e. those with no other pic- 
ture or boundary pixels neighbouring them, and all 
other pixels then have their charges set to 0. The re- 
sulting shower image is parameterised with a sec- 
ond moment analysis and for eventdisplay the time 
gradient across that image is found for the second 
pass of the charge integration. 

The image data from the individual telescopes is 
then tested against standard quality criteria for the 
number of tubes present in an image is > 5 (to en- 



sure a robust image axis), the minimum amount of 
charge in an image is size > 400 digital counts 
(providing run-to-run and NSB fluctation stability) 
and the angular distance in the camera is 0.05 < 
dist < 1.3 degrees (avoiding bias due to image 
truncation). If a sufficient number of telescope im- 
ages for an event pass the quality cuts the analysis 
proceeds to extract stereoscopic information like 
the direction on the sky and the shower core po- 
sition on the ground for each event. The image in- 
formation can then be compared to lookup tables of 
the parameters for simulated gamma rays as a func- 
tion of image size, impact distance to the shower 
core and the zenith angle of observation to achieve 
a greater rejection power of the uneven light distri- 
bution of cosmic ray events to that of the smooth 
compact gamma ray ones. Details of simulations 
of the VERITAS array can be found in [8]. For 
VEGAS the default method follows the prescrip- 
tion given in [3] where a mean scaled parameter 
(MSP) is given by 
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where corresponds to the parameter in ques- 
tion (width or length) for telescope i and 
Psim(0, size, r) is the mean value for simulations 
of a given size at a given impact distance (r) and 
zenith angle 8. Cuts based on the mean scaled pa- 
rameters for VEGAS are 0.05 < MSW < 1.02, 
0.05 < MSL < 1.15 (where W stands for width 
and L for length) and have a quality factor of 
Q r~j 1.8. In eventdisplay the normalized width 
method described in [7] is utilised to reduce the 
impact of outliers in the distributions. 
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where p sim (9, size,r) is the median value; and 
ago is the width of the distribution for 90% of the 
events. The cuts are -1.2 < MSCW < 0.5, 
-1.2 < MSCL < 0.5 with a quality factor of 
this form of cut is Q ~ 4. The values of the cuts 
given here should be considered standard for all re- 
sults presented, but are not necessarily considered 
as fully optimised. 
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Background estimation 

There is no set prescription for background esti- 
mation that can cope with every circumstance or 
source morphology. As such the best way to keep 
systematic uncertainties under control is by apply- 
ing several methods to the same data set and com- 
paring the results, a detailed description of com- 
monly used background estimation models can be 
found in [1]. 

Wobble mode observations [3] are favoured for 
point-like (or limited extension) sources since a 
background estimate can be derived from the 
events recorded during a run, allowing more on- 
source time and systematic effects in the back- 
ground estimation due to variation in weather or 
performance changes effectively cancel out. The 
background estimate is made from a region re- 
flected from the source position around the camera 
centre (the telescope pointing position). In order 
to gain a better estimate of the number of back- 
ground events multiple background regions can be 
used, each region set in a ring the same offset dis- 
tance from the camera centre, but still avoiding the 
region around the suspected source position to en- 
sure that poorly reconstructed gamma ray events 
do not bias the background estimate. Point sources 
show up as an excess of reconstructed shower di- 
rections close to the assumed source position. Fig- 
ure 1 shows just such a plot for 5 background re- 
gions on 4 good quality Crab nebula 0.5° wob- 
ble offset observations of 20 minute duration each. 
The spread of events closely matches that of a point 
source with an angular resolution of ~ 0.14° [8]. 
After a cut for a point like source of 9 2 < 0.025 for 
three telescope data (9 2 < 0.035 for two telescope 
data) the quality factor for all cuts is Q ~ 24. 

The reflected region methodology for wobble 
mode observations can essentially be applied to 
any part of the field of view displaced from the ob- 
servation position. This allows a 2-d map of events 
to be built up for the fov, figure 2 shows the 2-d dis- 
tribution of significances for the same Crab dataset 
using this method. A 2-d map like this will show 
up emission offset from the speculated source po- 
sition, sources of extended emission and a modest 
number of multiple sources in the same fov that a 
1-d analysis would miss. 
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Figure 1: 8 2 plot for 4 runs of Crab nebula ob- 
servations taken with a wobble offset of 0.5 de- 
grees is given by the line. The background esti- 
mate, made from 5 circular regions of 0.22 degree 
diameter each, is given by the crosses. 

An alternative 2-d mapping method is the ring 
background model. In this model a ring (in celes- 
tial co-ordinates) around a trial source position is 
used to give the background estimate. Since the 
ring covers areas with different offsets from that of 
the trial source position an acceptance correction 
function must be used in the normalisation for each 
position on the ring. Any part of a ring that crosses 
an assumed source position is also excluded from 
the background estimate. 

The source location accuracy can be seen in fig- 
ure 4 which shows the results for the fit of a 2-d 
Gaussian to the excess source counts. It can be 
seen that a source position can be accurately re- 
constructed to less than 0.05° for each run. 

Summary 

The standard VERITAS data analysis chain and the 
use of independent software packages to keep sys- 
tematic errors under control has been described. 
The use of different background estimation models 
allows hypothesis testing of different source mor- 
phologies and searching for unidentified sources 
within the field of view. The results from two dif- 
ferent analysis packages and three different back- 
ground estimation proceedures provide consistent 
results of ~ 30 a/y/hour for 3 telescopes on the 
Crab nebula for a wobble offset of 0.5°. 
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Figure 4: Reconstruction of the source position for a 
large number of Crab nebula runs (black squares). The 
red open cross is the overall reconstructed source posi- 
tion. 



Figure 2: 2-d sky map of significances with the re- 
flected region model for the same observations as fig- 
ure 1. 
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Figure 3: 2-d sky map of significances with the ring 
background model for the same observations as figure 1 . 
The white circle corresponds to the position of zeta tau. 
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